UNCLASSIFIED 


BOSTON  COLL  CHESTNUT  HILL  HA  DEPT  OF  PHYSICS  F/S  3, 

PROTON  EVENT  TIME  CHARACTERISTICS  AND  RADIO  BURST  DATA* (U) 

JUN  81  P  BAKSHX  t  T  N8UTEN  F19628-79-C-0084 

AF8L-TR-B1-0189  NL 


XAOO 


PROTON  EVENT  TIME  CHARACTERISTICS 
AND  RADIO  BURST  DATA 


AD  A  I  0 


rsi- 


Pradip  Bakshi 
Thao  Nguyen 


Boston  College 

Department  of  Physics 

Chestnut  Hill,  Massachusetts  02167 


June  1981 


Final  Report,  for  the  period  April  1979  to  February  1981. 


DTiC 


ELECTE 
NOV  2  1981 


A 


Approved  for  public  release;  distribution  unlimited. 


AIR  FORCE  GEOPHYSICS  LABORATORY 
AIR  FORCE  SYSTEMS  COMMAND 
UNITED  STATES  AIR  FORCE 
HANSCOM  AFG,  MASSACHUSETTS  01731 


®ll0300(5  7 


Qualified  requestors  may  obtain  additional  copies  from  the 
Defense  Technical  Information  Center.  All  others  should 
apply  to  the  National  Technical  Information  Service. 


SCCli'XiTy  CLASSIFICATION  OF  THIS  P»(,L  (  Hlifln  Dam  I  i.  lr..Ji 


/  REPORT  DOCUMENTATION  PAGE 


3ArP  READ  INSTRUCTIONS 

Mut- _ OEI-ORE  COMPLKTINf.  F'M« 

GOVT  ACL  C  5  SION  1:0.  ~  ftt.ClPLLN  T-S  CATALOG  NUMIJLIt' 


I.  HCPcrTT-SiH—BN  ■— -  |.\  GOVT  ACLCSSION  1:0.  J  fit  ClPU.  N  T  ‘5  C  AT  AlOG  NUMIJ  L II 

.  ATOLhTR-81-9169  [ _ - - 

«.  TITLE  (*b<t  SubtltlmJ  '  *  r  TYPE ’BI"WEV^«T  4  PERIOC  COVERED 

. .  - . .  /  r/  /  ~ 

Proton  Evtnt  Time  Characteristics  and  *..4  Final  /Repw't ' 

Radio  Burst  Data#  •  f  Aprtfl  1R79  Feba— ty  ^81 

•  .  —  >■  'MMtuwiin  mi 

7  AoTt-r. rc«;  b  'contract  on  gran’  num'ui'.T-^ 

Pradip  /feakshi  /  ,  «. 

Thao /Nguyen  j  /-■  F  19628-7^1-0084  ) 

T"~PER FORMING  ORGANIZATION  name  AMO  ADDRESS  '0.'  PROGRAM  ELEMENT  PROJECT.” 

AREA  6  WORK  UNIT  NUMBERS 

Department  of  Physics  62101  F 

Boston  College  /  /^  464303  AX 

Chestnut  Hill.  Ma  02167 _ >  /C  _ 

f».  CONTROLLING  DEVICE  N  AmE  ANU  ADDRESS  ^  ~  »2.  REPORT  DATE 

Air  Force  Geophysics  Laboratory  //  June  K81  / _ 

Hanscom  AFB,  Massachusetts  01731  —■  n‘  ufcn»twT»V'lp ag £1 

Monitor/William  R.  Barron/PHP _ 18 _ 

U.  MONITORING  AGENCY  NAtyfi^ft  AOCRESSfif  \htt*ronl  fr.*:n  i  nMrn/Jfnf  O/Rcr;  15  SECURITY  CLASS,  (of  thti  r»i>vti: 

/ *Y  /y  /  y  /  Unclassified 


«  .1  'WtRWqWliHHB  8 Mi  Rfcf>q«T-*t^  J E R^ 

b  'CON  TRACT  or  OMAN’  NUMUI'^l"; 

Ji  F  19628-7^-0084  / 

^  To”  PROGRAM  ELEMENT.  PROJECT  T  ASA " 
AREA  a  WORK  UNIT  NUMBERS 

62101  F 

<  /C  46-^03  A* 

L_-  ”  IJ.  REPORT  date  ” 

/}  ■  June  M81  / 

■—  la'  NbunrEwt)^  ag  is 
18 


^  ^//  Unclassified 


15a.  DECL  ASSlFlCATlON  DOWN'S  AOltiC 
SCHEDULE 


I  16.  DISTRIBUTION  STATEMENT  (ol  thin  Rtpnrt) 


Approved  for  public  release,  distribution  unlimited 


17  DISTRIBUTION  STATEMENT  ( ol  fh*  mb*tr*ct  entorod  in  Block  20,  II  dlllorent  from  Report 


in.  supplementary  notes 


KE  Y  WORDS  (C’.'ittnu*  on  reicrto  xidv  ft  nce*9*nty  mnd  i  l-ntlly  by  hh-ck  miTbrr) 

Solar  Protons  Predicting  Solar  Proton  Onset  Time 

Solar  Radio  Bursts  Predicting  Solar  Proton  Peak  Flux  Time 

Correlations,  Solar  Protons  vs  Radio  Bursts 


20  A9StRACt  (Cot'Urrj*  on  i»\-erso  iil>  If  nirci.t..?)1  O'xd  Mci;J»f\  b>-  t>f«»cfc  numberi 

^Various  studies  undertaken  for  relating  the  proton  event  time  character¬ 
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have  a  much  larger  onset  time  compared  to  the  Western  events.  For  western  hemi¬ 
spheric  events,  the  rise  time  of  the  radio  burst  is  found  to  be  well  correlated 
with  the  delay  before  onset  and  the  time  to  reach  the  peak  flux  level  of  the 
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I.  INTRODUCTION 


Solar  flare  radio  spectra  were  shown  to  be  predictors  of  significant 
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proton  events  over  a  decade  ago.  Hie  original  studies  provided  the 
basic  yes-no  criterion  (U  shaped  radio  spectrum).  Proton  event  magnitudes 
could  be  predicted  from  various  measures  of  the  radio  spectrum,  based  on 
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the  peak  or  mean  flux-density  and  the  rise  time  or  duration  of  the  burst. 

An  improved  technique  for  the  prediction  of  proton  magnitudes  has  also  been 
7-9 

developed,  based  on  (i)  using  the  radio  burst  energy  as  the  predictive 
measure  and  (ii)  applying  a  heliographic  longitudinal  correction  according  to 
the  flare  location.  The  radio  energy  parameter  is  obtained  by  a  time  integration 
of  the  radio  flux  at  a  given  frequency  for  the  duration  of  the  U-spectrum,  followed 
by  an  integration  over  a  broad  frequency  range.  The  spectral  character,  or  the 
slope,  of  the  proton  peak  flux  profile  was  shown  to  be  predictable  from  certain 
characteristics  (essentially  the  width)  of  the  radio  U-spectrum.10-1 2 

While  these  studies  enable  one  to  predict  the  proton  peak  flux  spectrum 
from  real  time  radio  burst  data,  they  do  not  provide  any  information  about 
the  onset  time  for  the  protons  or  the  time  when  the  proton  flux  will 
reach  the  peak  value,  or  more  generally,  about  their  over-all  time  development. 

The  present  project  was  initiated  to  tackle  this  broad  area,  and  in  particular, 
to  determine  some  radio  variables  or  characteristics  that  might  serve  as 
predictors  of  proton  time  characteristics.  This  Final  Report  describes  the 
various  efforts  undertaken  in  this  context. 

It  was  recognized  at  the  outset  that  proton  flux  data  with  good  time 
resolution  are  essential,  especially  for  the  study  of  onset  times.  The  IMP4 
proton  channel  29<E<94  MeV,  recorded  at  10  minute  intervals,  provided  the 
primary  data  base  for  this  purpose.  Proton  data  in  other  energy  channels 
were  also  compiled,  and  some  electron  data  were  also  recorded.  The  first 
study  dealt  with  the  effect  of  flare  location  on  the  onset  time  of  the 


particles.  As  described  in  section  II,  it  was  found  that  the  events  originating 
on  the  eastern  hemisphere  had  larger  time  intervals  between  the  occurrence  of 
the  flare  and  the  particle  arrival  time  compared  to  those  originating  on  the 
western  hemisphere.  However,  there  was  no  systematic  locational  difference  in 
the  onset  time  amongst  the  events  originating  on  the  western  hemisphere.  Thus 
the  data  base  consisting  of  only  the  western  hemispheric  events  would  provide 
the  proper  sample  for  the  study  of  correlations  between  radio-burst  characteris¬ 
tics  and  proton  flux  time  characteristics. 

Such  a  study  was  then  carried  out  using  the  rise  time  of  the  radio  burst 

(in  the  frequency  range  2-3  GHz  or  200-600  MHz)  or  the  rise  time  of  the  H^  flare 

as  one  of  the  correlation  variables  and  the  delay  before  onset  of  the  protons 

measured  from  the  peak  (or  start)  of  the  flare  or  the  radio  burst  as  the  other 

variable.  Good  correlations  were  obtained  as  described  in  detail  in  the 
13 

Scientific  Report  ,  and  as  summarized  in  Section  II.  It  was  also  shown  in 
that  report  that  the  delay  before  onset  of  the  protons  was  very  well  correlated 
with  the  delay  before  reaching  the  proton  peak  flux  and  also  that  the  latter 
time  interval  was  reasonably  well  correlated  with  the  radio  or  flare  rise  time. 
These  were  the  primary  results  obtained  in  this  program.  Numerous  secondary 
studies  were  carried  out  to  understand  the  remaining  "spread"  around  the  best 
fit  line  in  terms  of  other  parameters,  but  no  definitive  trends  were  found. 

The  main  conclusions  and  directions  for  future  work  are  discussed  in 


section  III. 
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II.  RESEARCH  WORK 


A.  Data  Base. 

The  primary  data  base  consisted  of  the  10  minute  interval  proton  flux 

data  in  the  29<E<94  MeV  proton  channel,  made  available  by  Smart  and  Shea 

in  digitized  and  semigraphical  computer  print  outs.  Typical  background 

-2  -1  -1 

levels  were  0.02  unit,  (unit  =  protons  cm  sec  ster  ).  The  actual 
onset  time  could  generally  be  determined  to  within  +  5  minutes.  The  delay 
before  onset  T^  is  defined  as  the  time  interval  between  the  peak  of  the 
radio  burst  (or  flare)  and  the  mid-point  of  the  ten  minute  time  interval 
during  which  the  proton  flux  rises  decisively  above  the  upper  range  of  the 
background  fluctuations.  A  related  quantity,  delay  before  onset  Tiq,  is 

lb 

defined  as  the  time  interval  between  the  start  of  the  radio  burst  (or  H 

a 

flare)  and  the  mid-point  of  the  ten  minute  time  interval  corresponding  to 

the  onset.  Similar  definitions  can  be  given  for  the  time  to  reach  the 

proton  peak  flux  level;  T„  if  measured  from  the  peak  and  T„„  if  measured 

M  Mb 

from  the  start  of  the  radio  burst  (or  H  flare). 

a 

All  the  events  of  interest  were  depicted  g-aphically  at  hourly  and 
10  minute  intervals.  The  peak  levels  ranged  from  0.1  to  several  hundred 
units.  Many  events  had  multiple  peaks  and  it  was  not  possible  to  define 
the  onset  time  unambiguously  for  some  of  them.  The  main  advantages  of  this 
channel  were  high  time  resolution  and  lack  of  electron  contamination.  The 
main  limitation  was  the  rather  restricted  time  range  (1967  to  early  1969) 
over  which  the  data  were  available. 

We  also  compiled  a  larger  data  base  (1967  to  1973)  using  other  proton 
channels  (>10,  >30,  >60  MeV)  and  electron  channels  (>40  KeV,  >2.7  MeV  etc.) 
from  the  computer  generated  comprehensive  catalog  of  particle  events  pre¬ 
pared  by  Smart  and  Shea.  Actual  onset  time,  peak  time  and  total  duration 


i 


I 
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in  each  channel  are  available  in  this  catalog.  Howe/er,  the  uncertainty 

in  these  times  is  much  greater  than  the +  5  minute  accuracy  of  the  directly 

read  data  in  the  29<E<94  MeV  channel. 

The  radio  and  flare  data  were  generally  taken  from  the  Simon  and 

Svestka^  catalog.  Other  sources^  ^  were  used  when  a  detailed  time 

history  or  other  information  was  available  for  some  of  the  events.  The  rise 

time  p  is  defined  as  the  time  interval  between  the  start  and  the  peak  of  the 

radio-burst  (or  the  H  flare) . 

a 

B.  Effect  of  Flare  Location. 

Each  hemisphere  was  divided  into  30°  zones  and  the  number  of  events  in 

each  zone,  n,  and  the  corresponding  average  delay  before  particle  onset  in 

hours,  T  ,  were  tabulated  for  each  of  the  years  1967,  1968  and  1969,  based  on 

14 

the  comprehensive  data  of  the  Simon  and  Svestka  catalog.  The  results  based 

on  all  events,  with  various  degrees  of  reliability  in  terms  of  the  flare- 

particle  association  are  displayed  in  Table  I.  If  we  restrict  the  data  to 

events  with  a  strong  flare-particle  association,  denoted  by  a  solid  circle  in  the 

Simon-Svestka  catalog,  we  obtain  Table  II  and  Figure  1.  It  is  clear  that  (a)  East 

zone  events  have  a  larger  delay  time  and  (b)  there  is  not  much  difference 

amongst  the  three  Western  zones  —  the  individual  scatter  of  events  is  greater 

than  the  differences  between  their  average  delay  times.  This  suggests  that  all 

West  zone  events  can  be  treated  as  a  group  to  study  what  other  parameters  are 

correlated  with  the  individual  delay  times.  It  should  be  noted  that  in  this  part 

of  the  study  the  particle  onset  was  not  defined  in  terms  of  a  particular  proton 

channel,  it  simply  refers  to  the  time  of  arrival  of  the  earliest  particles. 

We  also  note  that  the  results  obtained  here  are  in  broad  agreement  with  the 

18 

study  of  Reinhard  and  Wibberenz  . 


-5- 


TABLE  I 


1967 

1968 

1969 

Zone 

n 

T1 

n 

T1 

n 

T1 

E  90-60 

2 

5.0 

3 

5.84 

2 

8.75 

E  60-30 

2 

1.0 

10 

2.75 

2 

3.25 

E  30  0 

3 

1.0 

8 

3.06 

5 

4.30 

W  0-30 

8 

.74 

9 

.92 

5 

.41 

W  30-60 

11 

.54 

10 

.70 

8 

.45 

W  60-90 

16 

.60 

10 

.65 

13 

.44 

Total 

42 

TABLE  II 

50 

35 

1967 

1968 

1969 

Zone 

n 

T1 

n 

T1 

n 

T1 

E  90-60 

1 

2.0 

3 

5.84 

- 

- 

E  60-30 

1 

1.5 

5 

1.70 

1 

5.0 

E  30-0 

2 

1.0 

j 

2.84 

2 

3.5 

W  0-30 

4 

.68 

6 

.68 

5 

.41 

W  30-60 

3 

.35 

8 

.56 

6 

.53 

W  60-90 

3 

.58 

4 

.69 

8 

.47 

Total 

14 

29 

22 
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Correlations:  T.„  vs  p,  T,„  vs  T,„  and  Tu„  vs  p. 

_ lo _ Mb _ io _ Mb _ 

These  correlation  studies  are  the  main  subject  of  the  Scientific 
13 

Report  .  For  the  29<E<94  MeV  data  with  good  time  resolution,  the 

correlation  coefficient  for  T. _  vs  p  was  r  =  0.84  for  the  H  flare  and 

IS  a 

r  =  0.79  for  the  2-3  GHz  radio  burst  data  (Fig.  2).  Only  11  events  were 

available  due  to  the  limited  time  span  over  which  this  channel  was  operative 

and  due  to  our  selection  criteria.  The  results  for  other  energy  channels, 

but  using  the  same  events  as  above,  were  quite  similar.  When  the  data  base 

was  enlarged  for  the  other  channels  by  including  later  years,  the  correlation 

coefficient  was  somewhat  reduced  and  the  standard  deviation  somewhat  increased, 

but  the  best  fit  straight  line  remained  unaltered. 

For  T  vs  T,c  we  found  r  =  0.96  for  the  11  events  using  the  H  start 
MS  IS  a 

time  as  the  reference  time  and  r  =  0.94  using  the  2-3  GHz  start  time  as  the 
reference  time.  The  results  for  other  channels,  but  using  the  same  events, 
ranged  from  r  =  0.83  to  r  =  0.89.  With  the  data  base  enlarged  by  including 
later  years  the  correlation  coefficient  was  reduced  without  much  change  in 
the  best  fit  straight  line. 

For  T,„  vs  p  we  found  r  =  0.75  for  the  H  rise  time  for  a  data  base  of 
MS  a 

16  events  in  the  29<E<94  MeV  channel.  The  larger  sample  results  from  the 
fact  that  events  which  were  too  weak  tc  provide  a  reliable  onset  time  were 
quite  unambiguous  about  the  peak  time.  Similar  results  were  obtained  for 
the  2-3  GHz  radio  band  rise  time  correlations.  The  results  in  other  channels 
were  not  as  good,  due  to  their  poorer  time  resolution. 

Other  Studies. 

(1)  The  spread  around  the  best  fit  line  for  T^  vs  p  was  systematically 
examined  to  see  if  it  could  be  explained  in  terms  of  other  parameters.  It 


w.is  found  to  be  random  when  tested  against  location  (these  are  all  western 

hemisphere  events),  against  the  eventual  proton  flux  peak  magnitude, 

against  the  radio  spectral  shape  (U  or  non-U) ,  against  the  proton  spectral 

slope  and  against  the  width  of  the  U.  Similar  remarks  apply  to  T  „  vs  P 

MS 

spreads  as  well. 

(2)  The  peak  flux  level  I  and  the  ratio  R  =  I  /B  where  B  is  the  back- 

M  M 

ground  level  were  plotted  on  a  semilog  plot  vs.  T, ,  T,_,  T_,  and  1’  ,  but  no 

1  IS  M  MS 

clear  trend  could  be  discerned. 

(3)  Events  in  the  >40  KeV  electron  channel  were  examined  for  their 
onset  time  T^  vs  p.  There  were  many  more  electron  events  than  the  proton 
events.  If  we  considered  all  these  events  without  regard  to  their  peak  flux 
strengths  or  their  uncertainty  in  arrival  time,  the  sample  had  almost  35  events 
and  the  correlation  was  of  the  order  of  r  =  0.5.  In  most  of  these  events  the 
uncertainty  of  arrival  time  was  only  +  1  or  +  2  minutes. If  we  eliminate  the 
more  uncertain  (>  +  15  minutes)  events,  the  correlation  improves  considerably. 

If  we  further  eliminate  the  weak  events  (peak  electron  flux  <100  cm  ^  sec  1  ster 
the  data  base  is  reduced  but  the  correlation  improves  significantly.  Two  events 
v  'e  rather  special  and  there  seemed  to  be  problems  of  association  between  the 
flare  and  the  electron  increase.  These  have  been  excluded  from  the  above  con¬ 
siderations. 

(4)  Analyzing  the  proton  peak  flux  times  in  terms  of  a  two  step  process 

1 8 

represented  by  Tu  *  C,  +  C0/V  hws  been  suggested  by  Reinhard  and  Wibberenz  , 

M  1  Z 

where  is  a  longitude  dependent  propagation  time  and  C2  is  independent  of 
the  longitude.  The  second  term  represents  the  effect  of  velocity  dispersion. 

If  we  take  the  and  C2  for  the  various  western  hemispheric  events  as  given 


19 

in  Reinhard's  thesis  ,  we  can  test  whether  and  C2  individually  show 

any  correlations  with  p.  We  found  no  significant  correlations  between 

and  p  or  and  p ,  even  though  and  p  have  a  fairly  good  correlation. 

This  suggests  that  the  uncertainities  implicit  in  determining  (when  there 

is  not  enough  clear  velocity  dispersion)  bring  into  question  this  type  of 

separation  for  the  Western  hemispheric  events. 

(5)  The  decay  time  T  ,  defined  as  the  characteristic  time  in  the 

exponent  exp(-T/Tp)  for  representing  the  decay  phase  of  the  proton  flux  time 

profile,  was  noted  for  the  events  we  studied.  We  did  not  find  any  significant 

correlations  between  T  and  p,  T,_  or  Tw„ . 

D  IS  MS 

the  definition  of  T^  was  ambiguous. 


In  some  events  (with  multiple  peaks) 
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LII.  CONCLUDING  REMARKS 

The  main  conclusions  of  this  program  and  the  possible  directions  for  future 
work  are  as  follows. 

1)  There  are  good  correlations  between  T^s  and  p,  excellent  correlations 

between  and  T^,  ,  and  fairly  good  correlations  between  T^g  and  p.  These 

results  provide  a  prediction  scheme  for  T,_  and  T„_  in  terms  of  real  time  radio 

ib  Mb 

burst  data  for  western  hemisphere  events. 

2)  There  appears  to  be  a  single  characteristic  time  scale  for  each  event 
which  governs  the  radio  rise  time  the  proton  delay  before  onset,  the  proton  peak 
time  and  possibly  the  electron  onset  time  as  well.  Apart  from  the  "spread”,  one 
can  make  a  universal  model  for  the  time  development  of  the  various  components  of 
the  solar  flare  process;  empirical  observation  of  any  one  particular  component 
then  provider  the  time  characteristics  of  the  other  components.  In  particular, 

the  observed  (real  time)  radio  rise  time  can  serve  as  the  primary  factor  determining 
the  characteristic  time  scale  of  a  given  event. 

3)  Improvements  in  correlations  may  result  if  a  larger  data  base  were  available 
with  the  type  of  time  resolution  of  the  29<E<9A  MeV  channel.  Detailed  radio  burst 
time  profiles  at  various  frequencies  would  also  be  very  useful  in  judging  whether 

a  gi  radio  event  has  multiple  peaks,  and  if  so,  which  rise  time  is  the  appropriate 
one  to  use  for  the  correlations.  Improved  correlations  in  turn  would  lead  to  a 
better  prediction  scheme. 

A)  Further  work  with  various  electron  channels  would  be  valuable  in  judging 
which  channels  provide  high  correlations  with  the  radio  or  proton  time  scales.  If 
significant  correlations  can  be  shown  to  exist,  that  empirical  result  would  help 
guide  the  theoretical  modelling  of  the  flare  development  and  the  particle  propagation 
processes. 

5)  A  higher  radio  frequency  band  (e.g.  around  9  GHz)  can  be  used  for  the  radio 
rise  time  -stead  of  the  2-3  GHz  used  here. 
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6)  Besides  the  simple  rise  time  studied  here,  other  radio  parameters  such 
as  the  radio  peak  flux,  or  time  integrated  flux,  or  time  and  frequency  integrated 
flux,  the  width  of  the  U,  and  the  detailed  time  development  of  the  radio  spectral 
shape  should  be  carefully  examined  for  correlations  with  the  proton  time  scales, 

i 

or  for  understanding  the  spread  around  the  best  fit  lines  of  the  present  study. 

Our  preliminary  results  in  this  direction  were  negative,  but  a  larger  data  base 
is  needed  to  arrive  at  more  meaningful  results. 

7)  In  conclusion,  the  present  study  provides  a  scheme  for  assessing  the 
characteristic  time  scale  for  the  evolution  of  the  proton  event.  Earlier  studies 
gave  us  the  capability  to  assess  whether  a  proton  event  will  occur,  how  strong  it 
will  be  and  what  its  spectral  profile  would  be  like.  Here  we  have  taken  the  first 
step  to  assess  when  the  protons  will  start  to  arrive,  or  reach  the  peak  flux  level. 

From  the  practical  point  of  view,  knowing  the  time  scale  of  an  impending  event  is 
extremely  important. 
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Figure  1. 


Figure  2. 


FIGURE  CAPTIONS 


Particle  delay  before  onset  vs  Location  of  Flare 


Delay  before  onset  for  >29  MeV  protons  vs  Rise  time  in  the 
2-3  GHz  radio  frequency  band. 
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